Number of distinct sites visited by a subdiffusive random walker 
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The asymptotic mean number of distinct sites visited by a subdiffusive continuous time random 
walker in two dimensions seems not to have been explicitly calculated anywhere in the literature. 
This number has been calculated for other dimensions for only one specific asymptotic behavior of 
the waiting time distribution between steps. We present an explicit derivation for two cases in all 
integer dimensions so as to formally complete a tableaux of results. In this tableaux we include the 
dominant as well as subdominant contributions in all integer dimensions. Other quantities that can 
be calculated from the mean number of distinct sites visited are also discussed. 
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The mean number of distinct sites visited by a ran- 
dom walker on a lattice at time t after the start of the 
walk, <S'(t)iis a quantity often used to characterize such a 
walk ^^'^ many other quantities can be expressed 

in terms of it [3, §■ It is of course well known that S{t) 
depends on lattice dimension and geometry, and on the 
nature of the walk, that is, on the form of the waiting 
time distribution ^jj{t) between steps [3, II, 01 • ^ walk 
is said to be subdiffusive if the asymptotic mean square 
displacement of the walker grows sublinearly with time, 
e.g., (r'^it)) cx r with < 7 < 1, or (r'^it)) cx In'' t with 



/3 > [8|, l9| . One way to generate such subdiffusive walks 
is through waiting time distributions between steps that 
decay sufficiently slowly to possess no finite moments. In 
particular. 
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Here A is a dimensionless quantity and, in both cases, 
r is a constant with units of time and Y? is the mean 
square displacement of a single step. Noting the depen- 
dence of the mean square displacements on time, we will 
characterize the behavior of Eqs. ([T]) and ([2]) as slow and 
that of Eqs. ([3]) and ^ as ultraslow. 

The Laplace transform of will be denoted by 
(recognizable by its argument). For the inverse power 
law distribution one has [sj 
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and, for the logarithmic distribution, 
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In what follows we set r to unity. 

The probability Xn{t) that the walker has taken ex- 
actly n steps in time t is a multiple convolution over the 
^j{t) most easily expressed via the relation between their 
Laplace transforms [1, [ll| , 



LXn{t) = Xn{u)^[i,{u)T[l~ij{u)]lz 
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The mean number of distinct sites visited is given in 
terms of Xn{i) by 



S{r) = Y,SnXn{t), 



(8) 
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where S'„ is the mean number of distinct sites visited by 
the random walker in n steps. It then follows that 

00 

C [S{t)] = S{u) = [1 - i;{u)] ^" [^(")]" ■ (9) 

n=0 

The asymptotic behavior of S{t) for large t can be 
related via the discrete Tauberian theorem [3] to the be- 
havior of S{u) as M — > 0, that is, to the behavior of S{u) 
as ip{u) — > 1^. The discrete Tauberian theorem says that 
the expressions 
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are equivalent if p > 0, {c„} is a positive monotonic 
sequence, and F is slowly varying function at infinity 



2 



in the sense that F{\n) / F{n) ^ 1 as n — > cx) for each 
fixed positive A To evaluate S{u) as -u ^ and 
use this theorem to calculate the asymptotic behavior of 
S{t) we require the functional forms of Sn- These are 
well known [s!]: 
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For d = 2 and n — > oo 
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with the slowly- varying function at infinity, L[n), 
given by 
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and where 7 = 0.5792 ... is Euler's constant and 
a and b depend on lattice geometry. In particular, 
for a square lattice a — tt and 6 = 8. 

• For d > 3 and n — s- 00 the leading term is of 0{n) 
for all d but the subleading terms differ. 
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ay/n, d — 3 
alnn, d = 4 (15) 
a, d> 5, 



where the probability R of return to the origin 
depends on d, and where both R and the con- 
stant a depend on lattice geometry. In partic- 
ular, for a three-dimensional simple cubic lattice 
i? = 0.3405... 

Therefore, to evaluate S{u) it is necessary to evaluate 
sums of the form [see Eq. ^] 
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with u — > and f{n) ~ n", f(n) ~ Inn, and f{n) 
nL{n) for n 00. Defining 
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and applying the discrete Tauberian theorem one finds 
for u — > 0, 
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Using these results one finds that the asymptotic dom- 
inant and first subdominant contributions to S{u) for 
M — > are given by: 
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and 



■(aV^/2u)0i/2^ 



S{u) - (1 -^)^ + < (a/u)ln< 
a/u. 



d = 3 
d = A 
d>5. 



(23) 



To arrive at the time-dependent survival probabili- 
ties we need to apply the usual continuous Tauberian 
theorem 0, [l^ with a specific waiting time distribu- 
tion. For the inverse power law distribution ([T]) one has 
(j){l/u) ^ for M — > 0, so that the following asymptotic 
results for t — > cxD follow: 



S[t) 



72^/2 t-7/2 

r(l + 7/2) + V2r(l-7/2)' 



d=l 



(24) 



and 
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We have obtained these expressions via the straightfor- 
ward application of the procedure explained in [ll| com- 
plemented with the discrete Tauberian theorem. The 
first asymptotic term for d = 1 was obtained in [l^, [l4[ . 
On the other hand, for d = 2 we find the asymptotic 
(new) result 
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In particular, the dominant contribution, used in is 
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For the logarithmic waiting time distribution ([3]) one has 
(f){l/u) ^ Ahi'^{l/u) for M — > so that the new asymp- 
totic results follow: 



S{t) ~ V2lln^/2 t + In 
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where we have exhibited the third hnc to highhght the 
dominant term, and 

{(a%/l^/2) In'^/^ d = 3 
aln(^ln''i), d = 4 

a, d>5. 

(30) 

A number of observations about these results are in- 
teresting. First, we turn to the compactness of our ran- 
dom walks as measured by the ratio S{t)/V{t), where 
V(t) ^ (r^(t))'^^^ is the volume explored by the random 
walker. For the slow random walk this ratio behaves 
asymptotically as 
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For d = 1 the walk is thus compact (the walker visits 
every site in the region explored). In dimensions three 
or greater, the walk is non-compact. In d = 2 the loga- 
rithmic decay describes a walk that is "marginally non- 
compact." For the ultraslow walk we have 
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While in d = 1 the walk is again compact and in d = 2 it 
is again marginally non-compact, the decay marking non- 
compact behavior is slower here, and remains "marginal" 
(logarithmic) in all dimensions above d — 2. 



A second quantity of interest that can be calculated 
immediately from the number of distinct sites visited is 
the survival probability P(t) of a stationary target sur- 
rounded by a sea of moving particles of density p [ 
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The survival probability for the slow walk is a stretched 
exponential in all dimensions (with a logarithmic contri- 
bution in the exponent in two dimensions). Particularly 
interesting behavior is exhibited by the ultraslow walk. 
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^ln(ln'3 t), 
exp ^-cln'^ , d > 3, 

where c is a (different) constant in each case. Note that 
for d = 1 and (3 — 2 and for d > 3 and [3 — 1 the decay 
can thus be a power law, being even slower for smaller 
values of [3. As far as we can ascertain, this is a behavior 
not previously observed in the target problem. 

We have thus completed the panorama of "known" 
but previously not derived results (at least in the lit- 
erature that we could locate) for the asymptotic survival 
probability for an inverse power law waiting time dis- 
tribution ID) which leads to a slow growth of the mean 
square displacement. We have presented entirely new 
asymptotic survival probability results for the waiting 
time distribution ([3]), which leads to a growth of the 
mean square displacement characterized as ultraslow. We 
have concluded with a number of observations concerning 
the broader impact of these results, including the diverse 
range of slow decays of the survival probability in the tar- 
get problem associated with the ultraslow stepping time 
distribution. 
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